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Abstract 



Maxwell's first derivation of the equilibrium distribution function for a 

dilute gas is generalized in the spirit of the nonextensive g-statistics proposed 

by Tsallis. As an application, the g-Doppler broadening of spectral lines due 

to the random thermal motion of the radiating atoms is derived. 
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It is widely known that the thermodynamical or statistical description of nonextensive 
systems demand a generalization of the usual Boltzmann-Gibbs thermostatistics A 
few important examples of physical systems or processes where the standard approach seems 
to be inadequate are self-gravitating systems, some kinds of plasma turbulence, and self- 
organized criticality. 

The standard Boltzmann-Gibbs approach is based on the extensive entropy measure 

S = -k'^Pilnpi , (1) 

i 

where k is the Boltzmann constant and {pi} denotes the probabilities of the microscopic con- 
figurations. Ten years ago, in order to deal with the above mentioned difficulties associated 
with nonextensivity, Tsallis proposed the following nonextensive form of entropy 

where g is a parameter quantifying the degree of nonextensivity. For instance, given a 
composite system A + B, constituted by two subsystems A and B, which are independent 
in the sense of factorizability of the microstate probabilities (i.e. P/^^^^'* = Rj^^^Pj^^), the 
Tsallis measure verifies 

SqiA + B) = SqiA) + Sq{B) + (1 - q)S,{A)Sq{B) (3) 

In the limit g — > 1, 5*^ reduces to the standard logarithmic measure (1), and the usual 
additivity of entropy is recovered. There is a growing body of evidence suggesting that the 
g-entropy may provide a convenient frame for the thermostatistical analysis of many physical 
scenarios, such as stellar polytropes 0, turbulence in electronic plasmas 0, anomalous diffu- 
sion Levy distributions ll,!!^], the critical regime in low dimensional dissipative chaotic 
systems ||Tl]-|T^, the solar neutrino problem peculiar velocity distribution of galaxy 



clusters or more generally, systems endowed with long range interactions, long range 
memory effects, or a fractal-like space-time [jl],^. The nonextensive thermostatistics has 
been shown to be endowed with interesting mathematical properties ||TB|-0, the main the- 



orems of the standard statistics admiting suitable generalizations [|19| -|22|1. The issue related 



with the connection between g-statistics and g-thermo dynamics has also been addressed 
p.|23[], while the time evolution of Sg has been analyzed both in the discrete case through a 
direct application of the master equation and in the continuous one in connection with 



the Liouville and Fokker-Planck equations [25]. Some cosmological implications of Tsallis 



generalized thermostatistics have also been discussed pB| , p7| . However, at present there is 
only a limited understanding on the relation between the q parameter and the underlying 
microscopic dynamics. In the cases of low dimensional dissipative chaotic maps and 
in some toy models of self-organized criticallity , the value of q characterizing the system 



has been obtained from studies of the concomitant dynamics. In spite of the importance of 
these developments, they do not involve directly the Tsallis maximum entropy distribution 
and the experimental evidences supporting it. In order to clarify this point, let us briefly 
review the main observational facts supporting nowadays Tsallis' proposal. First, a Tsallis' 
maximum entropy distribution has been shown to describe properly a metastable equilibrium 
state of a 2-dimensional pure electron plasma 0]. Second, the g-distribution corresponding 
to an ideal classical gas provides a better fit for the observed distribution of peculiar veloci- 
ties of galaxy clusters than the ones obtained by recourse to more complicated models based 
on the standard thermostatistics [^. Finally, assuming a g- velocity distribution for the 
involved particles, the evaluation of the nuclear reaction rates in the solar interior predicts 
a neutrino flux in agreement with the observational data, thereby suggesting that Tsallis' 
thermostatistics may provide a solution for the well-known solar neutrino problem |]14| . 



It is remarkable that all the experimental evidence listed above deals, directly or in- 
directly, with the g-distribution of velocities, which can be obtained maximizing Sg under 
the normalization and mean energy constraints. Within a more general framework, such a 
distribution describes how the g-nonextensive canonical ensemble, associated with the clas- 



sical many body problem, depends on the particle velocities [^. In this way, the g- velocity 
distribution seems to be a reasonable nonextensive generalization of the celebrated Maxwell- 
Boltzmann distribution, which is recovered as the particular g — 1 limiting case. In spite of 
its theoretical interest, a satisfactory microscopic explanation for the physical origin of the 



g- velocity distribution is still lacking, although some interesting attempts have recently been 
made in connection with the linear and nonlinear Fokker-Planck equations [0. However, 
to shed some light on this matters, it seems important to consider suitable (nonextensive) 
generalizations of the kinetic approach pioneered by Maxwell and Boltzmann. 

In this letter, we are interested in exploring the kinetic route. Our aim is to rediscuss the 
correspondence between the parameter q introduced by Tsallis and the g-equilibrium velocity 
distribution for a Maxwellian gas, however, assuming from the very beginning a nonextensive 
generalization of the separability hyphothesis originally proposed by Maxwell . Hopefully, 



as happened in the extensive framework, this line of inquiry may provide some insight for 
a more rigorous kinetic irreversible treatment from the Boltzmann viewpoint. As a new 
application, we deduce a formula for the g-Doppler broadening of spectral lines. 

Let us now consider a spatially homogeneous gas, suposed in equilibrium at temperature 
T, in such a way that F{\')d^v is the number of particles with velocity in the volume element 
d^v around v. In Maxwell's derivation, the 3-dimensional distribution is factorized (lottery 
assumption) and depends only on the magnitude of the velocity pT|J5^ 



F (^^vl + vl + vi) d^v = f{v^)f{vy)f{v^)dv^dvydv^ , (4) 
from which it is straightforward to show that 

f{vi)=Aie ^ , i = x,y,z (5) 

and 

F{w) = A\e-^ , (6) 
where /5 = in order to recover the standard macroscopic thermodynamic relations, and 



Ai = y is the normalization constant. Naturally, in the nonextensive context described 
by (2), the starting basic hypothesis (4), which takes into account the isotropy of all velocity 
directions, must somewhat be modified. Physically, Maxwell's ansatz is tantamount to 
assume that the three components of the velocity are uncorrelated. However, this property 



does not hold in the systems endowed with long range interactions where Tsallis distribution 
has been observed []7|, p!^ , p!5| . Notice that the Maxwell ansatz is equivalent to express InF 



as the sum of the logarithms of the one dimensional distribution functions associated with 
each velocity component. A simple and natural way to generalize this procedure within the 
nonextensive formalism, in order to introduce correlations between the velocity components, 
would be to replace the logarithm function by a power law. However, in order to recover 
the ordinary logarithmic ansatz as a particular limiting case, it is convenient to express 
the power generalization in terms of the q'-log function, which is a combination of a power 
function plus appropriate constants. Elementary considerations may convince oneself that a 
consistent ^-generalization of (4) is (for simplicity we provisionally consider the bidimensional 
case) 



F d'v = e,{r\v.) In, /K) + r-\vy) In, f{vy))dv.dvy , (7) 

where the g-exp and g-log functions, eq{f), hig{f), are defined by 

e.(/) = [l + (l-g)/]^/^-^ , (8) 

K/ = ^^Y^ • (9) 

As one may check, eg(lng /) = lng(eq(/)) = /, and from (9) we see that the g-log differentia- 
tion, ^ lug / = f~^^ is also satisfied. Note also that in the limit q ^ 1 the identities (8)- (9) 
reproduce the usual properties of the exponential and logarithm functions, and (7) reduces 
to the bidimensional case of (4), as should be expected. Now, partial g-log differentiation of 
(7) with respect to Vx yields 

d\n,F _ d\n,[e,{p-\vx)\n, fjvx) + p-\vy)\n, fjvy))] 
or equivalently, 



(10) 



where x = y'^x + "^y ^ prime means total derivative. Note that analogous equations 
apply to the remaining components even whether we had considered the n-dimensional case. 
Introducing the shorthand notation 



1 F'ix) 



(12) 



we may rewrite (11) as 

Hx) 



l__d_ 

Vy dVy 



(13) 



Vx dvx 

The second member of the above equation only depends on v^, while the third one is a 
function exclusively of variable Vy. Hence, equation (|13]) can be satisfied only if all its 
members are equal to one and the same constant, not depending on any of the velocity 
components. So, we may put $(x) = —mj, where m is the mass of the particles and 7 is 
an arbitrary constant. Of course, the introduction of m at this point is dictated only by the 
known Maxwellian limit. As one may see from (9), /^~^(^^x) l^q /(f^) = lug* f{v^), where 
q* = 2 — q. Hence, the general solutions for f{vx) is given by (equivalent expressions are 
valid for f{vy) and fivz)) 



lUg* f{Vx) 



(14) 



where, without loss of generality, we have written the integration constant in a convenient 
form. Now, taking the q-exponential in both sides of (14) it follows that 

,2\i 1/1-9* 



f{Vx) 



l + {l-q*)\ hv A 



(15) 



and defining a new constant /3 as 



7 



7 



1 + (1 - q*) lUg. A A^-1* 
we find the generalized expression 



(16) 



f{vx) = Ag 



l-(g-l)^ 



1 

q-l 



(17) 
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where we have introduced a subindex q to make exphcit the g-dependence of A. From 
(17) we see that the Gaussian probabihty curve of the Maxwelhan gas is replaced by the 
charactheristic power law behavior of Tsallis' nonextensive framework, and as expected, the 
limit 9 = 1 recovers the exponential extensive result. Note also that for values of q greater 
than unity, the positiviness of power argument means that (17) exhibits a thermal cut-off in 
the maximal allowed velocities. The components of the velocities lie on the interval [— L], 



where L = y ^^^^^i^ ■ Hence, the integration limits in the standard normalization condition 
is modified in such a way that only if g = 1 they go to infinity. Taking this into account one 
may show that the normalization constant Ag can be expressed in terms of Gamma-functions 



as 



^^"l^J r(^) ■ ^^^^ 

Further, using that limi^i^oo ^-^^e"'^^"^ = 1 (see Abramowitz [^]), it is easy to see that 
Ai is the standard Maxwellian result. 

Before continuing we need to obtain the complete distribution. By adding the Vz com- 
ponent to (7) it is readily seen that 

F {^^vl + vl + vfj d'^v = eg(ln,. f{v^) + In^. f{vy) + In^. f{vz))dv^dvydvz . (19) 

Hence, taking the q-logarithim of the above expression and repeating the same algebraic 
steps of the one-dimensional case it follows that 



q-l 



(20) 



where Bg is fixed by the 3-dimensional normalization condition. We find 



3 

m 



2-KkT 



(21) 



As expected, the q-distribution (20) is isotropic meaning that all velocity directions are also 
equivalent in this generalized context. As remarked earlier for the one- dimensional case, 
there also exist a temperature dependent cut-off on the magnitude of the velocities. From 



(18) we see that the Bi — Af — [2^]^ is the standard 3-dimensional MaxweUian result as 
it should be. 

Now we discuss an important point of principle, which is related with the definition of 
marginal probabilities in the context of g-statistics. In ordinary space, a marginal probability, 
say, f{vx), may be obtained from the 3-dimensional distribution by 

f{vx) ^ j F{v^,Vy,v^)dVydv^ . (22) 

This elementary and natural definition, widely applied in statistical physics, is not usually 
regarded as deserving any further scrutiny. However, within the nonextensive framework 
we are discussing here, the concept of marginal probability distributions shows some new 
remarkable features. Basically, this occurs because the distribution f{vx), as given by (22), 
does not coincid with our equation (17). As one may check, it has a power different from 
In other words, using the above formula, the 3-dimensional distribution (20) leads to a f{vx) 
with a different value of q. Strictly speaking, since the q-distribution cannot be factorized, 
the power of the one-dimensional distribution obtained from a marginal probability like (22), 
depends on the number of spatial dimensions. The power increases by | for each additional 
dimension present in the complete distribution. This is equivalent to say that f{vx) is also a 
power law, but with an effective g-parameter. However, to interpret this fact in a consistent 
way, it is required to introduce an effective temperature in the one-dimensional distribution. 
Naturally, the same happens when we consider arbitrary dimensions m, n, where m < n. 
In particular, this means that the zeroth law of thermodynamics is not satisfied for systems 
described in this nonextensive framework. Hopefully, a proper modification of (22) can 
be found which avoids these undesirable features, but preserves the interesting ones. In 
this concern, we suggest the following general expression for the marginal m-dimensional 
distribution in a n-dimensional g-velocity space 

J,, X J F''{vi,V2,...,Vn)dv„,+ ldVm+2- ■ -dVn 

f{Vi,V2,...,Vm) ^ -— — — , (23) 

J F°'{Vi,V2, ...,Vn)dVidV2...dVn 

where a = 1 — ^{q — l){n — m) . For q = 1, we have a = 1, and the standard definition is 
recovered. In the above discussed case, n = 3, m = 1, one finds a — 2 — q. For this value 

8 



of a it is straightforward to show that (23) reproduces the one-dimensional distribution 
(17). As a matter of fact, our prescription (23) solves the conflict with the zeroth law of 
thermodynamics, in the sense that the same power law and temperature of the complete 
distribution is always obtained regardless of the specific dimensions involved in the problem. 
Besides, it is interesting to realize that equation (23) can be interpreted as defining the 



ordinary marginal probability function computed using the escort distribution ||34 



(24) 



instead of being evaluated using the original distribution F . This is a standard procedure 
in the fractal thermodynamic formalism |34 . 



Broadening of Spectral Lines. The random motion of particles broadens spectral lines, 
first, because of collisions between the particles (pressure broadening), and second, due to 
the thermal Doppler effect of the radiating atoms. As widely known, in the extensive case 
the first effect is proportional to pT~^, where p is the pressure, whereas the second scales 
with Ta. Let us now discuss the latter effect using the g-Maxwellian velocity distribution. 



The standard result was derived by Lord Rayleigh and further observed by Michelson ||35| , |36 

In order to estimate the magnitude of the g-Doppler broadening for the visible light 
emitted by the molecules of a hot gas, it is enough to consider the one dimensional case. 
Neglecting relativistic effects, the frequency shift viewed along the x direction is given by 
the standard Newtonian formula 

z/ = 1^0 (l + ^) • (25) 

The frequency distribution expected for a spectral line centered at vq is obtained changing 
variables from Vx to v. From (17) it is readily seen that 

9-1 g 



1-iq-l) 



-dv , (26) 



2kT V 

which has also the form of a Maxwellian g-distribution. The broadening is usually measured 
by the width of the spectral line at half intensity. It is easy to check that in terms of the 
wavelenght, the standard deviation is replaced by 
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AA^ = An 



2kT (I -2^''^' 



1/2 / , ^ N 1/2 

-Aj^21n,2 . (27) 



For g = 1 the standard result is recovered as should be expected p2| , p5| . However, although 
mantaining the same temperature dependence, the thermal Doppler broadening is mod- 
ified in the nonextensive framework. Naturally, the above formula can be used to limit the 
q-parameter. Note that in a log-log plot of AA^ versus the straight line is displaced 
parallel to itself for each value of q different from unity. In comparison with the standard 
g = 1 result, the above velocity distribution gives a narrowing of the Doppler width for 
g < 1, and a broadening tendency for q larger than unity. 

Tsallis' Generalized Mean Values. It is worth noticing that we might have considered 
a different generalization for the factorization condition. In principle, instead of equation 
(0) one may assume the following ansatz 

F = e,{\n, f{v,) + In, f{vy)). (28) 

In this case, repeating the same steps we have already explained, instead of (^\), one obtains 
a slightly different velocity distribution 

1 



1- (1 



1-0 

(29) 



which is exactly that one determined by Tsallis MaxEnt prescription when the generalized 
mean value 

(vl), = I fvldv, (30) 



is a meaningful constraint. On the other hand, by employing the standard linear mean 
values, distribution ( [T7| ) is obtained. As a matter of fact, none of the main results and 
conclusions of this paper change in a significative way whether one adopts the alternative 
factorization prescription (|28|). For instance, all the results corresponding to our analysis of 
the g-Doppler broadening follow simply by replacing everywhere {q — 1) by (1 — g). This 
is a strong indication that a more conclusive result needs a full kinetic theoretic treatment, 
which requires a proper generalization of the Boltzmann if-theorem. 
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It is possible that the approach developed here may be implemented even if more general 
expressions for nonextensive entropies, which reduce in a commom limit to the Boltzmann- 
Gibbs-Shanon form, are considered. Tsallis measure Sq is not the only conceivable mathe- 
matical generalization of the standard logarithmic entropy. Other interesting nonextensive 
entropic functionals have been recently proposed [^,^. However, Tsallis entropy has been 
shown to be endowed with many elegant and useful mathematical properties, and its associ- 
ated g-MaxEnt distributions have been experimentally observed [0,|l5l . It would be of great 
interest to explore the, so far poorly known, mathematical properties of the recently intro- 
duced entropies ||37| , |38[| , as well as to determine if they admit relevant physical applications. 



Conclusions. In the present work we have obtained Tsallis non-extensive velocity distri- 
bution by recourse to an argument akin to the celebrated derivation advanced by Maxwell 
for the equilibrium velocity distribution. As shown by Maxwell, his distribution is the only 
one compatible with isotropy and factorizability with respect to each velocity component. 
Similarly, Tsallis q-distribution is uniquely determined by the requirements of (i) isotropy, 
and (ii) a suitable generalization of the factorizability condition. Maxwell's factorization 
condition is tantamount in requiring that the logarithm of the complete distribution func- 
tion be equal to a sum of terms, each one depending only on one velocity component. 
Instead of logarithm, our factorization condition requires that a power of the complete dis- 
tribution be equal to a sum of terms, which also depend only on one velocity component. 
Reformulating this last condition in terms of the Tsallis g-logarithm function. Maxwell's 
expressions are recovered in the g — 1 limit. The same happens with the formula giving the 
broadening of spectral lines. 

It is important to stress that the simple transformation (g — 1) — (1 — g) is enough 
to recast all the present results within the complete Tsallis formalism |^ based on Tsallis 
generalized mean values. In that case, Tsallis g-distribution would adopt the form 



F(v) = B, 



1-9 

(31) 



where Bq is given by the expression obtained from (plD after replacing (g — 1) by (1 — g). 
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In a similar vein, the g-Doppler width of spectral lines would appear under the guise 



kT 



1/2 



AArf = AJ^2''lng2 . (32) 

Mathematically, Maxwell factorization condition is similar to the Maxwell-Boltzmann 
hypothesis of "molecular chaos" , which assumes that the two-molecule distribution function 
describing the colliding molecules is factorizable as the product of two one-molecule distri- 
butions (i.e. F(vi,V2) = /(vi) /(V2)). This hypothesis plays a fundamental role in the 
standard kinetic theory of gases . Boltzmann' proof that any initial velocity distribution 



evolves irreversibly towards Maxwell's distribution does not rely just in the general princi- 
ples of classical mechanics. It also needs the additional assumption of "molecular chaos". 



Boltzmann himself |^9[ recognized that the hypothesis of "molecular chaos" may not always 
hold, especially at high densities. Our present results suggest that Boltzmann's approach 
to Maxwell's velocity distribution can be adapted to the non-extensive setting by recourse 
to an appropriate generalization of the "molecular chaos" assumption. This issue will be 
addressed in a forthcoming communication. 
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